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ABSTRACT 

To date the search for habitable Earth-hke planets has primarily focused on nuclear burning stars. 
I propose that this search should be expanded to cool white dwarf stars that have expended their 
nuclear fuel. I define the continuously habitable zone of white dwarfs, and show that it extends from 
«0.005 to 0.02 AU for white dwarfs with masses from 0.4 to 0.9 Mq, temperatures less than « 10* 
K, and habitable durations of at least 3 Gyr. As they are similar in size to Earth, white dwarfs 
may be deeply eclipsed by terrestrial planets that orbit edge-on, which can easily be detected with 
ground-based telescopes. If planets can migrate inward or reform near white dwarfs, I show that 
a global robotic telescope network could carry out a transit survey of nearby white dwarfs placing 
interesting constraints on the presence of habitable Earths. If planets were detected, I show that the 
survey would favor detection of planets similar to Earth: similar size, temperature, rotation period, 
and host star temperatures similar to the Sun. The Large Synoptic Survey Telescope could place even 
tighter constraints on the frequency of habitable Earths around white dwarfs. The confirmation and 
characterization of these planets might be carried out with large ground and space telescopes. 
Subject headings: astrobiology — binaries:eclipsing — eclipses — planetary systems — planets and 
satellites:detection — white dwarfs 



1. INTRODUCTION 

The search for habitable planets has focused on stars 
similar to the Sun as it is the sole example we have of 
a star with a habitable planet, and nuclear burning pro- 
vides a long-lived source of energy (jKasting et al.lTl993l : 
iLunine et al1l2008[ ). White dwarfs, which are as common 
as Sun-like stars, may also provide a source of energy for 
planets for gigayear (Gyr) du rations. Whi te dwarfs have 
typical masses 0.4—0.9 Mq (jProvencal et al. 1998), but 
have radii onl y ft!l% of the Sun, about the same size 
as the Earth ( Hansenll2004f) . The most common white 
dwarfs have surface temperatures of ~5000 K which are 
referred to as "cool white dwarfs" since hotter white 
dwarfs are easier to detect (iHansen 2004). Cool white 
dwarfs typically have luminosities of 10~'^Lq, so a planet 
must orbit at «0.01 AU to be at a temperature for liq- 
uid water to exist on the sur face, the so-called habit- 
able zone (jKasting et al.l[l993[ ). The small size of white 
dwarfs can cause large transit depths by Earth-sized 
or even smaller bodies which could in principle be de- 
tectable with ground-based telescopes (iDi Stefano et al.l 
I2OIOI: iFaedi etaIll20lol: iDrake et al.ll2blOD . 

Prior to becoming a white dwarf, a Sun-like star ex- 
pand s to a red giant, eng ulfing planets within «1 AU 
(e.g.. lNordhaus et al.ll201 O). Planets present in the white 
dwarf habitable zone (WDHZ) must arriv e after this 
phase . This may occur via several paths (|Faedi et al.l 
l2010f ) : planets can form out of gas near the white dwarf, 
via th e interaction or merger of binary stars (jLivio et al.l 
l2005fl. or b y capture or migration from larger distances 
(jPebes fc~S igurdsson 2002). There are precedents for 
each of these processe s: one neutron star has a plane- 
tary system (Wolszczan fc Fraill 119921) which may have 
been formed from a disk created after the supernova 
jPh innev fc HansenI 119931 ): pulsars show low-mass stel- 
lar companions being whittled down to planet masses 



(|Fruchter et al.|[l988() : and white dwarfs show infrared 
emission and atmospheric compositi ons indicative of 
close orbiting dust and/or bodies (e.g.. lZuckerman et al.l 
120101 ). Consequently, short period planets around white 
dwarfs might plausibly exist. I sidestep the question of 
formation, which has had little theoretical attention, and 
instead address the location and duration of habitable 
zones around white dwarfs (Section [5]), the detection of 
planets in the WDHZ, even if their frequency were much 
less than 1% (Sections|3]and|4l), and how characterization 
of these planets might proceed (Section [5|) . 

2. WHITE DWARF HABITABLE ZONE 

I co mpute the WDHZ boundary following the proce- 
dure in lSelsis et al.l (|2007[ ). I determine the white dwarf 
luminosity and effective temper ature, T^s, versus ag e 
from white dwarf cooling tracks (jBergeron et al.l [20011) . 
With the luminosity and effective temperature, I com- 
pute the range of distances within which an Earth-like 
planet could have liquid water on the surface if it were 
placed there with an intact atmosphere. I use compu- 
tations of the limits of the habitable zone for stars of 
different effective temperatures that are based on em- 
pirical limits from our solar system combined with one- 
dimensional radiative — convective atmospheric models 
for Earth-like planets that include water loss at the in- 
ner edge an d the maximum CQ2 greenhouse effect at the 
outer edge (jKasting et al.llT993[) . 

Most white dwarfs have masses Mwd — O-QMq and 
CO interior composition, for which I plot the WDHZ ver- 
sus time in Figure 1 as a blue shaded region. This region 
shrinks with time as the star cools. A planet enters at the 
bottom of Figure 1 and moves vertically up the figure as 
its white dwarf host ages, so it starts off too hot for liquid 
water, passes through the WDHZ, and then becomes too 
cold. The duration a planet spends within the WDHZ, 
tuz, has a maximum of 8 Gyr at «0.01 AU. Based on the 
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Fig. 1. — WDHZ for Mwd = O.GMq vs. white dwarf age and 
planet orbital distance. Blue region denotes the WDHZ. Dashed 
line is Roche limit for Earth-density planets. Planets to right of 
dotted line are in the WDHZ for less than 3 Gyr. Planet orbital pe- 
riod is indicated on the top axis; and white dwarf effective temper- 
ature on the right axis. Luminosity of the white dwarf at different 
ages are indicated on right. 
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Fig. 2. — CHZ vs. white dwarf mass and planet orbital distance. 
Green region is the CHZ for tjjz > tmin=3 Gyr, H-atmosphere. 
Left solid line is Roche limit for Earth-density planets. The other 
lines show how the CHZ outer boundary changes for tniin=l Gyr 
(dotted line), for tmin=5 Gyr (dash-dotted line), or for an He at- 
mosphere with imin= 3 Gyr (dashed line). Horizontal line indicates 
the most common white dwarf mass of 0.6 Mq, plotted in Figure 
1. 



WDHZ limits, I next define the "continuously habitable 
zone" (CHZ) as the range of planet orbital distances, a, 
that are habitable for a minimum duration, imin (Figure 
2). I choose a minimum duration of tmin = 3 Gyr, which 
results in a CHZ within a < 0.02 AU: a planet that or- 
bits within this distance spends at least 3 Gyr within the 
WDHZ. From 0.4 to 0.9 Mq with tniin=3 Gyr the outer 
boundary of the CHZ always falls within 10% of 0.02 
AU for hydrogen and helium atmospheres (Figure 2). To 
check the sensitivity to the white dwarf cooling compu- 
tations, I have also compu ted the WDHZ with BASTI 
models (jSalaris et al.l[2010[ ). for which I find a slightly 
longer ^hz- I have also computed the CHZ for atmo- 
spheres with t,„in = 1 Gyr and 5 Gyr which shifts the 
outer boundary of the CHZ outward/inward by a factor 
of Ril.5/0.7 (Figure 2). 



There are several important consequences of the 
WDHZ and CHZ. First, the range of white dwarf tem- 
peratures in the portion of the CHZ within the WDHZ 
is that of cool white dwarfs, «3000-9000 K (right hand 
axis in Figure 1), similar to the Sun. At the hotter end 
higher ultraviolet flux might affect the retention of an 
atmosphere, these planets would need to form a sec- 
ondary atmosphere, as occurred on Earth. Excluding 
higher temperature white dwarfs only slightly modifies 
the CHZ since they spend little time at high temper- 
ature. Cool white d warfs are photometrically stable 
([Ftantaine fc Brassarll l2008l ) , which is critical for find- 
ing planets around them. Second, for white dwarfs 
with temperatures >4500 K, the WDHZ is exterior 
to the Roche limit (tidal disruption radius), or = 
O.OO54AU(pp/pe)-i/3(MwD/O.6Af0)i/3, where pp,® are 
the mass densities of the planet and of Earth. Conse- 
quently, the 3 Gyr CHZ lies between or < a < 0.02 AU, 
indicated with the green region in Figure 2. 

Third, the CHZ occurs at white dwarf luminosities 
of 10~^'^ to 1O~^L0, about 10 magnitudes fainter than 
the Sun, which sets the minimum telescope size for de- 
tection. Finally, the orbital period of white dwarfs in 
the CHZ is w4-32 hr. At this period the timescales 
for tidal circularization and tidal locking are wlO-lOOO 
years , so rocky planets w ill be synchronized and circular- 
ized (jHeller et al.ll201ll ): the side of the planet near the 
star will have a permanent day, while the far side will 
have a permanent night. The planet will orbit stably as 
it cannot raise a tide on the compact white dwarf. Radi- 
ation drag will not cause the orbit to decay, but magnetic 
field dr ag might, depending on the planet's conductivity 
(iLTeT al. 1998). 

3. DETECTION OF PLANETS IN THE WHITE DWARF 
HABITABLE ZONE 

I plot example light curves of planets in the WDHZ in 
Figure 3 for sizes within a factor of « 2 of Earth orbiting a 
white dwarf near the peak of the white dwarf luminosity 
function. The transits last w2 minutes, and have a max- 
imum depth of 10% — 100%; these events can be detected 
at a 100 pc distance with aim ground-based telescope. 
To detect planets in the CHZ one must monitor a sample 
of white dwarfs for the duration of the orbital period at 
0.02 AU, and planets present in edge-on orbits will be 
seen to transit their host stars. The transit probability 
is 

Rp/R^ + i?wD/O.O13i?0 



Ptrans 



1.0% 



a/O.OlAU 



(1) 



so for every 100 Earths orbiting white dwarfs at 0.01 AU 
with random orientations, on average one will be seen to 
transit. 

I define ry© to be the number of planets with O.IM® < 
Mp < lOMe in the 3 Gyr CHZ (a < 0.02 AU). To mea- 
sure to an accuracy of Ri33% requires detecting w9 
planets, so one must survey ^lO'^rj^^ white dwarfs. 

4. SURVEY STRATEGY 

The local density of wh ite dwarfs is (4.7 ±0.5) x 10~^ 
pc"-^ ([Harris et al.l 120061) . For a survey out to Dmax < 
200 pc, the white dwarfs should be nearly isotropically 
spaced on the sky at one per w2(100pc/r'max)^ deg^. 
This exceeds the field of view of most telescopes, so each 
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Fig. 3. — Example light curves of habitable Earth-like planets 
transiting a 0.6 Mq white dwarf (T^ff = 5200K, _RwD = O.OlSiJ©, 
inclination = 89°. 9, 10% linear limb darkening, and a = 0.013 
AU) with masses of O.lMg (red, about Mars mass), l.OM^ (black, 
Earth twin), and lO.OAffg (blue, super-Earth). The ratios of the 
planet radii, Rp, to white dwarf radius, -RwDi indicated. 

white dwarf must be separately surveyed for the presence 
of transiting planets. 

I have simulated an all-sky survey with a worldwide 
network of 1 m aperture telescopes to monitor the white 
dwarf CHZ (typically 32 hr, during which telescopes 
distributed in longitude follow a single star) following 
[Nutzman & Charbonncau (2008) to compute the tele- 
scope sensitivity, including sky and read noise, and as- 
suming an exposure time of 15 s; I conservatively ex- 
panded the error bars an additional 50%. Each white 
dwarf is given a multi-planet system whose innermost 
planet is chosen from a log-normal centered at 2aR with 
width of 0.5 dex, with subsequent planets packed as 
closely as allowed by dynam ical stability on a timescale 
of 10^ yr (jZhou et al.ll2007D . leading to a uniform distri- 
bution in log a, with a gradual cutoff within 2aR. The 
planet masses are drawn from dn/dMp cx M~°' from 

10-2 Afe (wMoon) to lO^M® («Saturn), with a = 4/3 
to match the observed slope measured with the Kepler 
satellite (|Borucki et al.ll201ll ). as well as that found in 
the solar system. I assumed Earth-composition planets 
with the rad i us det ermined from the mass according to 
iSeager et al.] (pOOl . Two exposures within transit with 
a signal-to-noise of at least 6 constitute a detection; this 
keeps the false-positive level to <1% for the entire sur- 
vey. I then scale the simulated detection rates with 77^ 
to determine the expected number of detected planets. 

To detect 9±3 planets, for rj^ = 50% a survey of 2800 
white dwarfs within 52 pc is required for a total on-sky 
time of 10 years (for a single telescope with 31 hours 
per white dwarf on average) . A smaller planet frequency 
of = 10% requires «20,000 white dwarfs within 100 
pc for a total 69 years of telescope time on sky. For a 
network of twenty 1 m telescopes distributed around the 
globe, such as th e Las Cumbres Ob servatory Global Tele- 
scope (LCOGT: iHidas et al.l[200l. or the Whole Earth 
Telescope (WET lNather et al.lll99Ql ), devoted to observ- 
ing white dwarfs at 25% efficiency (50% of time at night, 
50% weather loss), the total calendar time required would 
be 2 years for a survey of 2800 white dwarfs and 14 years 
for a survey of 20,000. If the CHZ is surveyed to only 



0.01 AU, this would decrease the required calendar time 
to 8.5 months and 5 years, respectively, but would also 
decrease the planet yield. 

Figure 4 shows the probability distribution for plan- 
ets detected in 10"* survey simulations. Each one sur- 
veys 20,000 single white dwarfs out to 100 pc for planets 
within 0.02 AU. Of the detected planets, an average of 
40% will be currently within the WDHZ. Remarkably, 
the detection probability peaks near planets of the size 
and temperature of Earth due to the coincidence in size 
of the Earth and white dwarfs, and the coincidence be- 
tween the WDHZ at the peak of the white dwarf luminos- 
ity function and 2aR. This leads to the following biases: 
(1) large planets have a higher transit detection proba- 
bility oc Rpdn/ dRp, so the number detected declines if 
dn/dRp is steeper than R~^; (2) small planets cause shal- 
lower transits for Rp <C i?wD ~ -Re > so the detection rate 
scales as oc R^dn/dRp (jPepper et al.ll2003l ). although the 
range of luminosities of white dwarfs flattens this decline; 
(3) cooler planets have a smaller probability of transit, 
so fewer are detected as oc Tp for small Tp and a uniform 
distribution in log a; and (4) hotter planets orbit hotter 
stars, whic h are less num erous, dn/dTwD Twd^ for 
large Two (|Hansenll2004[ ). Although these trends should 
occur for any volume-limited survey, the break in the ra- 
dius detection limit depends on the size of the telescope 
and signal-to-noisc cutoff that is chosen: larger telescopes 
or smaller signal-to-noise cuts will be more sensitive to 
small radius planets. Figure 4 is sensitive to the prop- 
erties of the planet population: if the inner cutoff, Oin, 
is further /closer, then the peak temperature moves to 

cooler/hotter temperatures, oc Cj^^ , and the total num- 
ber of planets detected declines/increases as a'^^ due to 
the lower/higher transit probability. If the planet size 
distribution is steeper/flatter, then the detected planet 
distribution peaks at smaller/larger sizes. 

Prior to such a survey, a nearby sample of cool white 
dwarfs must be found using measurements of the re- 
duced proper motion. Ongoing and planned deep as- 
trometric surveys, such as 90 Prime, Skymapper, Pan- 
STARRS, URAT, and GAIA, should find most cool 
white dwarfs out to 100 pc (V<2\) w ithin the decade 
(jHenrv et al.l 120091: iKalirai et all |2009[) . Some of these 
surveys might also find transits if the requirement of 
three epochs in transit is relaxed. For example, GAIA 
will observe 200,000 di sk white dwarfs 50-100 times each 
(jPerrvman et al.l 120011) . possibly detecting one epoch in 
eclipse for «10% of these stars with habitable transiting 
planets. 

For values of 770<1O%, more stars must be observed to 
detect planets, so a better strategy is to observe multi- 
ple white dwarfs simultaneously with a wide-field imager 
with fast readout, suc h as the L arge Synoptic Survey 
Telescope (LSST LSST Science B ook 2009). The LSST 
survey is expected to detect wlO'' white dwarfs over half 
of the sky at >5<t to r<24.5 with wlOOO epochs each 
and two 15 s exposures per epoch over the duration of 
the 10 year survey. Since LSST is a magnitude-limited 
survey, the white dwarf temperature distribution peaks 
at lO'^K. I have taken the si mulated dete cted distribu- 
tion of white dwarfs for LSST (jJuric et al.ll2008.) , created 
simulated light curves for white dwarfs with planets, and 
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Fig. 4. — Probability density, d^?i/((iTpdlog_Rp), of detected 
planets vs. planet radius, Rp (log axis scale), and planet effec- 
tive temperature, Tp (assuming the same albedo as Earth). The 
contours enclose 25%, 50%, and 75% of all detected planets; the 
contour levels are 29%, 53%, and 76% of the peak density. Earth 
(ffi). Mars (d"), and Venus (5) symbols indicate the radii and effec- 
tive temperatures of these solar system planets. 



added noise (|LSST Science Bookl[2009l) . I find LSST can 
detect >9 GHZ planets if r]Q>5xlO~^, where detection 
requires that at least three epochs fall within transit with 
two points each detected at >7cr. The LSST survey will 
be biased toward detecting shorter period (oc P~'^/^) and 
large-size planets that have yet to enter the WDHZ since 
their stars are hotter. This could be improved by either 
continuously observing some fields for several nights, or 
taking more exposures per field, resulting in detection 
of smaller, cooler planets, thus constraining smaller ry^. 
LSST will identify the white dwarfs with reduced proper 
motion measurements as the survey is being carried out. 

I estimate that >10^ double white dwarf eclipsing bi- 
naries with orbital periods similar to WDHZ planets will 
be found with LSST us ing the BSE popul ation synthesis 
model for binary stars (jHurlev et al."2002') with parame- 
ters taken from observed binaries (Ragliavan ct al. 2010). 
I find that about 2.5% of white dwarfs will have a white 
dwarf companion with a period in the range of 8 — 64 hr 
which might be mistaken for a transiting GHZ planet if 
these are viewed edge-on. Follow-up of planet candidates 
will be required to distinguish the two possibilities: white 
dwarf binaries will show primary and secondary eclipses 
of different depths (if the two white dwarfs differ in tem- 
perature), offset secondary eclipses due to light travel 
time, gravitational lensing (jAgoll |2002[ ). and Doppler 
modulation, and an eclipse shape that differs from plan- 
etary transits if non-grazing. I simulated light curves of 
white dwarf binaries including these effects and find that 
in the worst case of two white dwarfs with identical tem- 
peratures these distinguishing features may be detected 
with 10 — 100 m ground-based telescopes for systems out 
to 50 — 100 pc, either photometrically or with radial ve- 
locities. Another concern is grazing eclipses from white 



dwarf/M dwarf binaries; these can be identified by the 
eclipse shapes, spectral energy distribution, and differing 
secondary eclipse depths. 

5. PLANET CHARACTERIZATION 

The parallax and spectrum of a white dwarf yield 
its mass, luminosity, atmospheric composition, and ra- 
dius; then the transit depth gives the planet radius. 
The planet's mass cannot be measured from Doppler 
shifts due to the featureless spectra of cool white dwarfs, 
but may be bracketed by the range of com positions for 
planets of a given size ISeager et al.l 120071) . The mass 
might be measured by observing wavelength dependent 
absorption, such as Rayleigh scattering that varies as 
8H/Rp In A ()Lecavelier Pes Etangs et al.ll2008[ ). where H 
is the atmospheric scale height, causing the transit depth 
to vary by wfew millimagnitude over 400 — 500 nm if 
Rayleigh scattering dominates over other sources of opac- 
ity. If atmospheric mol ecular weight can be es timated, 
so can the planet mass ()Miller-Ricci et al.ll2009[ ). If two 
planets transit a white dwarf (about 2.5% of the time for 
mutual inclinations within 5° and in a packed planet sys- 
tem), then transit timing variations might constrain the 
planet mass es if the orbital period ratio is nearly corn - 
mensurate (jHolman fc Murray! [2001 lAgol et al.l 120051) . 
Infrared phase variation of planets ( Knutson et al.ll2007l) 
in the GHZ is R:i0.5%— 2% at 15—20 ^m; however, I es- 
timate that the James Webb Space Telescope {JWST) 
cannot detect this for cool white dwarfs due to telescope 
noise. For a hotter white dwarf of 10^ K at 100 pc with 
an Earth-like planet at 0.01 AU with a day — night con- 
trast of 30%, the phase variation of 0.1% at 7.7 ^m might 
be detectable at 9a with the MIRI imager on JWST 
(I Swinvard et al.ll2004D . 

Several topics require further study. The global cli- 
mate models for determining the WDHZ should include 
fast synchronous rotation, magnetic fields, varied planet 
atmospheric composition, radiogenic heating, and white 
dwarf cooling. The WDHZ is a necessary but insufficient 
criteria for habitability. For example, planets that start 
hot may not retain their atmospheres, as has also bee n 
argued for planets orbiting M dwarfs (jLissaueil I2007D : 
this may require vo latile delivery from more distant bod- 
ies in the system (jJura &: Xull2010D or planetary out- 
gassing. To retain an atmosphere might require a larger 
planet escape velocity, possibly favoring super-Earths for 
habitability. 

Formation mechanisms must be modeled to help mo- 
tivate future surveys. For example, gravitational in- 
teractions of a planet and star with a third compan- 
ion body m ay be responsible for creating hot Jupiters 
(iFabrvck y" & Trcmainc 2007), which is also promising for 
moving distant planets around white dw arfs to 2aR, w 
0.01 AU, the tidal circularization radius (jFord fc Rasid 
|2006[) . It is also possible that tidal disruption of 
a planet or a companion star will result in the for- 
mation of a disk whic h may cool and form planets 
(jGuillochon et all I2010D , out of which a second genei- 
atioii of planets might fo rm (jMenou et al. 2001; Pereti 
[20T(il: IHansen et al.l[20()9t ) . 

The most common white dwarf has T^s ~5000 K, close 
to that of the Sun; consequently, inhabitants of a planet 
in the GHZ will see their star as a similar angular size 
and color as we see our Sun. The orbital and spin pe- 
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riod of planets in the CHZ are similar to a day, caus- 
ing Coriolis and thermal forces similar to Earth. The 
night sides of these planets will be warmed by advec- 
tion of heat from their day sides if a cold-trap is avoided 
(jMerlis &: Schne ider 2010j. Transit probabilities of hab- 
itable planets are similar for cool white dwarfs and Sun- 
like stars, but the white dwarf planets can be found us- 
ing ground-based telescopes (e.g., LCOGT, WET, and 
LSST) at a much less expensive price than space-based 
planet-survey telescopes. 
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Note added in p roof. Rene Heller informed me that 
iMonteirol (|2010 f) discussed the white dwarf habitable 
zone for known white dwarf stars, showing a boundary 
similar to that in Figure 1. 
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